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Abstract -At the FERMILAB TEVATRON and CERN LHC accelerators partonic processes 
reach such a high energy region, where the weak interactions could enter into consideration in 
the hadron yields through the opening of resonant W and Z exchange diagrams. We investigate 
these contributions to inclusive charged hadron production in p + p and p + p reactions at a few 
I I TeV energies, and compare our results to existing CDF and CMS data. We perform a leading 

order perturbative QCD calculation and include electroweak channels to estimate the order of 
magnitude of charged hadron excess. The energy dependence of hadron production and the fine 
^ ^ structure of transverse momentum distributions are investigated. As expected, in p + p collisions 

^ electroweak channels yield much larger contribution, than in p+p collisions. However, these types 

<^ of contributions remains close to be negligible at the available accelerator energies. 
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Introduction. — The first year of successful run of 
CERN Large Hadron Collider resulted enormous amount 
of data \n p + p collisions at energies Vs = 2.36, 7 TeV. 
The preliminary data of the ALICE ^, ATLAS [2] and 
CMS [3] collaborations on inclusive charged hadron pro- 
duction display the extension of transverse momentum 
spectra well measured up to 100 GeV/c or even beyond. 
These results complement results from the CDF Collabo- 
ration at FERMILAB TEVATRON. They reported (4)[5] 
the inclusive charged hadron distribution measured up to 
PT « 50 - 100 GeV/c at ^ = 1.96 TeV collision energy. 
With increasing luminosity and running time we can ex- 
pect sufficiently large number of recorded events at LHC 
so that the inclusive hadron spectra will be measured up to 
transverse momenta of couple of hundred GeV with high 
precision. Thus all contributions beyond standard chan- 
nels based on strong interaction particle production will 
be revealed, if they have small but non-negligible effect. 

Inclusive hadron production at high transverse mo- 
menta is used to be calculated in the framework of pertur- 
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bative quantum chromodynamics (QCD) improved parton 
models, which reproduces the well known scaling proper- 
ties of the measured hadron spectra, namely E^^s^ oc Py", 
where n w 5. When electroweak channels from W and Z 
boson exchange start to give noticeable contributions then 
one would expect the above scaling to be violated because 
of the resonance structure of the heavy boson exchange. 
This violation appears at half of the bosonic mass values 
at jet levels, but shifted down to lower momentum values 
in hadronic spectra as a consequence of jet fragmentation. 

In this letter we address a general question. The moti- 
vation of our work is to quantify the role of electroweak 
channels in inclusive hadron production. Recently, the 
CDF collaboration [4| reported a very strong scaling vi- 
olation in inclusive charge hadron spectra at high trans- 
verse momentum, which triggered considerable theoretical 
interest. In Ref. [6] the fragmentation functions have been 
investigated carefully, and no means have been found to 
explain these CDF data. In another effort the inclusion of 
electroweak channels were proposed to explain the CDF 
data [7]. Here we report the outcome of our attempt to 
repeat this calculation. Although we make several approx- 
imations, we carefully investigate their possible effect and 
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Fig. 1: The dominant weak-interaction mediated Feynman di- 
agrams with a hadronic final state. 



conclude that both at TEVATRON and at LHC the yield of 
the electroweak channels remains below the per cent level 
as compared to the yield of the strong processes. 

In the meantime, the CDF coUaboration pubhshed a pa- 
per with reanalyzed data on inclusive hadron spectra |5], 
invalidating their earlier findings. The strong scaling vi- 
olation reported before almost entirely disappeared, al- 
though the new data stop at smaller transverse momentum 
value (pt < 50 GeV/c), than the firstly reported dataset 
{PT < 100 GeV/c). Also recently, the CMS collaboration 
has published their data at ^/s — 7 TeV and no strong 
violation was reported. Thus, it seems that strong scaling 
violations do not appear in the inclusive hadron spectra at 
current colliders. Nevertheless, we believe that for future 
reference it is useful to publish our quantitative predic- 
tions with the electroweak production channels also taken 
into account. 

Pure perturbative QCD framework. The in- 
clusive hadron spectra is determined in the perturbative 
QCD improved parton model 



X fa/p{Xa,Q'^)fb/p{p)(xb,Q'^)Dh/c{Zc,QF) (1) 
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where all reactions with u, d, s or c quarks in the incoming 
or outgoing channels are summed over. Here fa/p denotes 
the parton distributions function of parton a in the pro- 
ton, and Dji/c is the fragmentation function of hadron h 
emerging from parton c. For factorization scales we use 
Q = KpTjct = kpt/zc, and Qp = npx, with k = 1 as 
defauh. We use the MSTW2008 parton distribution func- 
tions [9] and the KKP fragmentation functions 10 for our 
leading order (LO) calculation. 

Estimating the vireak-interaction contributions. 

Here we estimate the order of magnitude of the weak- 
interaction contributions by using a leading order calcu- 
lation in the parton model. The dominant W boson ex- 
change channels are QaQb — >■ qcqd, as illustrated by 
Fig. [ij We neglect Z boson exchange and i-channel dia- 
grams. This is justified by the fact that the coupling of Z 



to quarks is smaller than the W coupling, and t-channel 
diagrams are suppressed as compared to the resonant s- 
channel contributions. We also neglect photons, as these 
do not posses a resonance structure, therefore their con- 
tribution is also expected to drop rapidly with px- 

The Born-level cross section for the partonic subprocess 
shown in Fig. [l]is 
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where s and u are the partonic Mandelstam variables, mw 
and Tw are the mass and width of the W boson, V is 
the CKM matrix, and — e^ /sin^Ow is the weak cou- 
pling. Since the main contribution comes from s values 
close to the W mass, we can use the weak coupling mea- 
sured at the vector boson mass, neglecting the running 
of the electroweak coupling, which has weaker momentum 
dependence, than a^. Thus we use a{Mz) = 1/127.9, and 
the weak mixing angle sin^Owi^z) = 0.231 (see 11 ). 



Difference between pp and pp collisions. Con- 
sidering usual strong interaction production channels the 
hadron yields are very similar in pp and pp collisions. The 
slight deviation in the parton distribution function for va- 
lence and sea quarks will generate a very small difference 
in the production of charged hadrons. 

In earlier investigation of on-shell and W~ produc- 
tion, an asymmetry has been found |12|, which strongly 
depends on the d ~ u asymmetry in proton distribution 
functions. Similar effect appears in weak boson mediated 
parton-parton scattering. We expect a larger contribution 
from -channels in pp collisions relative to the yield in 
pp collisions because in the p the weight of valence anti- 
quark is much larger than the the weight of the sea anti- 
quarks in p. 

Numerical results. — We calculate hadron produc- 
tion for p+p colhsions at ^/s = 1.96 TeV, and compare the 
numerical results to the CDF data (4j[5]. We repeat the 
investigation in our pQCD based calculation at 7 TeV in 
p + p collisions and compare the numerical results to CMS 
data. The factorization scale dependence was investigated 
by varying the k values between 1/2 and 3/2. Figure 2 and 
3 display the obtained results. The data points are from 
CMS f3l, and CDF [4j|5] respectively, which are close to 
the pQCD LO predictions. 

Our results indicate that weak channel contribution is 
more than three orders of magnitude smaller then strong 
channel contributions. The choice of the scales do not 
influence significantly this result. 

One would expect a bump in the jet spectrum at px ~ 
from the decay of W bosons, and a bump in the 
hadron spectrum at somewhat lower momenta, due to 
fragmentation of quarks to lower momentum hadrons. In- 
deed our calculation of the weak-interaction contribution 
reproduces this feature. 
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Fig. 2: [Color online.) Contributions to inclusive average 
cliarged particle yield from strong (upper red band) and weak 
(bottom blue band) interaction channels in pp collisions at 
■^s — 7 TeV energy. The scale uncertainty of the strong and 
weak interactions is indicated by the bands. Experimental data 
are from the CMS Collaboration |3j. 



Fig. 3: (Color online.) Contributions to inclusive charged par- 
ticle yield from strong (upper red band) and weak (bottom blue 
band) interaction channels in pp collisions at -y/s = 1.96 TeV 
energy. The scale uncertainty of the strong and weak interac- 
tions is indicated by the bands. Experimental data are from 
the CDF Collaboration [5j. 



As mentioned earlier, the weak-interaction contribution 
should be considerably higher in p + p collisions than in 
p + p. To investigate this difference, we made calculations 
at different center of mass energies. We define the ratio: 
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and plot this quantity on Fig|4] During this study, we 
fix K = 1. The main difference in the curves for different 
a/s energy values is that in higher energy collisions the 
ratio Rpp/pp is smaller. This can be readily understood by 
the increasingly larger weight of the gluon in the proton 
with decreasing partonic momentum fraction x because 
the higher the collision energy, the smaller x contribute 
to the produced hadron with the same transverse momen- 
tum. The ratio increases with increasing px due to the 
increasing valence quark contribution, coming from the 
higher momentum fraction part of the parton distribution 
functions. 

We also define the ratio of the weak and strong interac- 
tion contributions in pp collisions: 

Rw/s = \ li"^, : -• (4) 



We investigate the dependence of this ratio on the col- 
lision energy and the transverse momentum. Our result 
are displayed on Fig. 5. We also used k = 1 here. As 
expected, the ratios have a peak at somewhat lower trans- 
verse momenta than Mw/2. The ratio is higher in lower 



energy collisions because of the dominance of gluons at 
small X, discussed above. 

We study the sensitivity of the above results on the 
scale choice k, and find that varying k between 1/2 and 
3/2 makes the ratios vary a few ten percents (which is 
the usual sensitivity of leading order QCD results on the 
factorization scale) . The main conclusion, namely the neg- 
ligibility of weak interactions, remain the same. 

Conclusion. — We have investigated the inclusive 
hadron production in pp and pp collisions at FERMILAB 
TEVATRON and CERN LHC, including strong and weak- 
interaction channels. We have performed a LO calculation 
to obtain an estimate on the hadron excess from weak 
channels. 

Our parton model calculation confirms the usual prac- 
tice of neglecting weak interaction contributions when cal- 
culating hadron spectra, in pp and pp collisions at col- 
lision energies accessible with present technology. The 
weak-interaction contribution is always orders of magni- 
tude smaller. 

We expect that repeating our analysis at next-to-leading 
order accuracy would lead to similar conclusion especially 
for the ratios Rpp/pp and Rw/s where the effect of the 
NLO K factors cancel to a large extent. 

A naive estimate of comparing the running couplings 
at the momentum scales involved, namely {asw /cts)'^ = 
O(10~^) is still an overestimate of the ratio on Fig. 5, be- 
cause of the dominance of gluon channels in hadron pro- 
duction. 
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Fig. 4: {Color online.) The ratio of weak-interaction contribu- 
tions in p -|- p to p -|- p collisions at different -/s energy values. 
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